EFFECTS OF H 2 0 VAPOR ON VIBRATIONAF 
REFAXATION IN EXPANDING AND CONTRACTING FLOWS 


Willard E. Meador,* Lawrence W. Townsend' and Gilda A. Miner* 
NASA Langley Research Center 
Hampton, VA 23681-0001 


Abstract 

As opposed to previous explanations based 
on the effects of anharmonicity of simple 
diatomic molecules, traces of water vapor are 
suggested to be the most likely cause of the 
anomalously fast vibrational relaxation of such 
gases observed in supersonic and hypersonic 
nozzles. The mechanism is the strong V-VR 
coupling with H?0 molecules that dramatically 
facilitates the collisional transfer of vibrational 
energy. Slight moisture content is thus a real- 
world aspect of gas dynamics that must be 
considered in characterizations of shock tubes, 
reflected shock tunnels, and expansion tubes. 


Introduction 

A well-known problem in nonequilibrium 
diatomic gases is the greatly accelerated 
vibrational relaxation observed in expanding 
supersonic flows compared with theoretical 
predictions based on relaxation times deduced 

from shock-tube data. 1 ' 6 Since this discrepency 
is impossible to explain with molecules modeled 
as linear harmonic oscillators (LHO's), the 
anharmonicity of real molecules is usually 
invoked to provide the rationale. 
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The qualitative description of the effects of 
anharmonicity is very simple. As originally 
proposed by Treanor, Rich, and Rehm, 1 a gas 
expands through a sequence of vibrational 
population distributions given by the product of 
the Boltzmann ratio of number densities in the 
r th and ground-state levels of LHO molecules 
and the Treanor correction factor 
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where f>,- is a characteristic temperature defined 
such that k A ; . is the energy of the r th level and T 
is the translational temperature. Equation (1) is 
most easily derived by assuming that rapid, near- 
resonant V-V energy exchange processes enable 
corresponding detailed balancing on a time scale 
short compared with that of the drop in T. 
Accelerated relaxation in an expanding gas is 
then seen to be a two-step phenomenon: (1) 

V-V pumping of upper levels because the 
anharmonic difference (r Tf | - f>,-) in Eq. (1) 
increases with increasing r; (2) the operative step 
of enhanced V-T relaxation rate coefficients for 
molecules in these upper levels because of the 
narrower spacings between levels. The starting 
or stagnation temperature of the expanding flow 
must be at least comparable to fli for the 
Boltzmann factor not to kill the participation of 
upper levels. None of the conditions for 
accelerated relaxation is met in contracting 
gases, as caused, for example, by the passage of 
a shock wave. 

As compelling as the qualitative arguments 
for anharmonicity effects may seem, 
quantitatively they can account for only a small 
fraction of the accelerated relaxation observed in 
the expanding flows in reflected shock tunnels. 



Calculations by Rich and Treanor, 3 and more 

recently by Ruffin and Park, 4 predict 
enhancements by at most a factor of two, 
whereas typical data show relaxation occurring 
up to two orders of magnitude faster than 
allowed by LHO molecules. The experimental 
stagnation temperatures and departures from 
equilibrium are not large enough for 
anharmonicity to do more. 

Clearly, then, some additional effect is 
needed to explain the relaxation data. The 
purpose of this paper is to suggest that traces of 
water vapor in an otherwise pure diatomic gas 
may be enough to account for the differences. 


Binary Mixtures 

A fundamental part of any description of 
nonequilibrium processes is proper treatment of 
the tendency to approach equilibrium demanded 
by the second law of thermodynamics. In 
particular, a partial differential equation is 
needed to describe the behavior of the average 
vibrational energy e* per molecule after a 
system has been knocked out of equilibrium, 
either by compression or by expansion. 
Compression deposits energy initially into the 
translational mode, leaving vibrational levels 
underpopulated relative to full equipartition of 
energy; expansion and cooling leave the 
vibrational levels relatively overpopulated. 
Subsequent behavior can be described 
mathematically by implementing the methods of 
reaction kinetics. 

For present puiposes, we assume no 
chemical reactions (including dissociation), no 
multiple quantum jumps in a molecule due to 
inelastic collisions, and no contributions from 
conventional transport processes. We also 
assume LFIO molecules since it is effects other 
than anharmonicity that we are now seeking. 
Straightforward application of rate equations for 
the populations of individual vibrational levels 
yields the following expression for the average 
vibrational energy per molecule of species i: 



where the subscript e denotes equilibrium values 
at the same translational temperature (introduced 


through detailed balancing) and the x's are 
relaxation times. 

The same application of rate equations that 
gives Eq. (2) also gives the relaxation times in 
terms of rate coefficients k ( (defined herein with 

the units atmosphere' 1 x second" 1 ). To be more 
specific, we now consider the binary mixture 
N 2 -H 2 0. Let subscripts 1 and 2 refer to N 2 and 

FI 2 0, respectively, when applied to e* , X, 
characteristic vibrational temperature f>, and 
mole fraction x. Relevant reactions are 


N 9 + N 2 — (k : )— > N 2 + N 2 ; (3a) 

N 9 + H 2 0 — (k 2 j— > N 9 + H 2 0 ; (3b) 

N 9 + H 2 0 — (k 3 )— > N 9 + H 2 0 ; (3c) 

H 9 0 + H 2 0 — (k^) — 4 H 2 0 + H 2 0; (3d) 
H 2 0* + N 2 — (k 5 )-» H 2 0 + N 2 . (3e) 


The rationale for the critical reaction (3c) is 
given in the next section. 

An important feature of Eq. (2) is the last 
term, which couples the relaxation of species i 
with those of the other species and depends on 
the existence of exchange reactions like (3c) 
above. Its physical interpretation is clear: it 

represents impedance of the relaxation of 
species i by the relaxation of other species. If 
e\ < e\ e and e* 2 <e* 2e , for example, species 2 

can increase its vibrational energy not only by 
V-T collisions with all molecules, but also by 
V-V collisions with molecules of species 1 that 
leave the latter with less vibrational energy and 
thus headed in the direction opposite from 
where they need to go. This particular V-V 
scenario and its reverse will contribute to each 
of Xp X 12 , t 2 , and x 2 i in Eq. (2). 

A point to be made here is that CFD codes 
treating all reactions as V-T would misrepresent 

the essence of N 2 * deexcitation. In particular - , 
(3c) would be missing from the above list and 
the entire burden of deexcitation of N 2 * by 
collisions with H 2 0 would fall on the V-T 
reaction (3b). This would mean two things: 
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(1) a totally nonphysical value of k 2 that may 
have to be several orders of magnitude too large 
in order to compensate for the neglect of ky, (2) 
a mathematical format without the explicit 
coupling between Eq. (2) for N2* and Eq. (2) 

for H2CE. Applications of such codes would be 
questionable. 

The principal relaxation times X,- in Eq. (2) 
can be written as follows for N2-H2O with 
reactions (3a)-(3e): 
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Special Properties of H2O Molecules 

Very much unlike the homonuclear N2, 
water molecules have strong dipole moments to 
cause rapid V-T relaxation by reaction (3e), 
involving both dipole-quadrupole and dipole- 
induced dipole interactions, and especially by 
the self-deexcitation of (3d). Accordingly, and 
no matter what the relative mole fractions, the 
comparatively small X2 in Eq. (4b) enables the 
average excitation energy of PbO to achieve a 
quasi-steady state (QSS) at each step of the 
much slower N2 relaxation. More quantitatively, 
the QSS's must exist if the inequality 
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is satisfied, which appears to be true for all N2- 
PbO mixtures and also for all Cb-tEO mixtures. 
Good experimental evidence for the QSS’s is 
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provided by acoustical absoiption in moist N2. 7 

Implementation of the QSS's means 
De 2 / Dt = 0 and thus 
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where k 3 = k 3 exp 
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is the rate 


coefficient for the backward V-V reaction. 


Apart from the k 3 and k 3 terms, these results are 
typical of combination rules usually employed 
in CFD codes. The strictly coupling relaxation 
times are 
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from Eq. (2) for PQO. Elimination of \e 2e -e* 2 
from Eq. (2) for N2 then yields 
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Dt x 



( 8 ) 


in which the subscript 1 has been dropped and 
the effective single relaxation time for N2 
satisfies 


1 

x 


X 1 X 12 X 21 


(9) 


1 

P X 21 


X 1^3^2 


( 

1 

V 



\-l 


'\,e 


k$ 


•(5b) 


1/ 


Recognition of the QSS's thereby allows the 
mathematical formulation of the N2-H2O binary 
mixture to be reduced to that of a single-species 
gas. 
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Water molecules have another special 
property that is even more important for present 
purposes. Their unique geometry of two light 
hydrogen atoms bonded roughly at a right angle 
between them to a far more massive oxygen 
atom means an unusually small moment of 
inertia. This, in turn, means easily induced 
rotational transitions and very closely spaced 
rotational energy levels; consequently, excellent 
energy matches abound between vibrational 
levels of most molecules and vibrational- 
rotational levels of H2O. That water vapor, even 
in very small amounts, is observed to be a superb 
quencher of vibrational excitation should not be 
surprising. 

Under most laboratory conditions, the rele- 
vant levels for the exchange reaction (3c) are the 
first excited vibrational state ofN2 ($1 = 3395 
K) and the corresponding portion of the near- 
continuum of rotational levels superimposed on 
the first excited state of the bending mode of 
FhO (f>2 = 2295 K). Fairly crude collision 
modeling of this reaction, which is more aptly 
termed near-resonant V-VR instead of obviously 
nonresonant V-V, was performed by Nagel and 

Rogovin. 8 These authors make the key point 
that conventional selection rules do not apply to 
hard collisions, which means that significant 
jumps in rotational quantum numbers can and 
do occur in N2-H2O exchange collisions. They 
also achieved reasonable agreement with 
experimental data. 

The bottom line is that the rate coefficient k 3 
for reaction (3c) is not negligible; in particular, 
it is certainly orders of magnitude larger than k \ , 
which represents the only reaction occurring in a 
pure N2 gas where vibrational relaxation is 
extremely slow because of the noninteractive 
nature of homonuclear molecules (no dipole 
moment). Accordingly, the comparatively few 

N 2 -H 9 0 collisions in a slightly moist N2 gas 

will be far more effective in quenching N2 

* 

vibrational excitation that the many N 2 -N 2 
collisions. Such is the overwhelming signifi- 
cance of traces of water vapor in understanding 
nonequilibrium phenomena in laboratory test 
facilities. 


Correction Factors for Slightly Moist 
Diatomic Gases 

From this point on, we assume the following 
sequence of rate coefficients obtained from 


acoustical measurements 7 of moist N2 and air 
and corresponding to reactions (3a)-(3e): 


> > k 5 > k 3 > k 3 > > k 1 and k 2 . (10) 

The respective values at room temperature are 
10 9 , 1.2 x 10 6 , 1.2 x 10 5 , 3.2 x 10 3 , and 0.082 
(atm.-sec.) _1 . 

Though not used specifically in the brief 
present analysis, we also note that the rate 
coefficient for 


o 2 + h 2 o — > o 2 + h 2 o* 


is about two orders of magnitude larger than the 
£3 in Eq. (3c) because of the close match 
between fHC^) = 2239 K and fKFbO) = 2295 
K. Flere the near-resonance actually is V-V 
instead of V-VR. 

Substitution into Eq. (9) of Eqs. (4a), (4b), 
(5a), and (5b), together with use of the 
inequalities in Eq. (10), gives the following 
approximate expression for the effective N2 
relaxation time appropriate for Eq. (8): 


P x P*i { k& ij 
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( 11 ) 


The more revealing notation h = X 2 to denote 
humidity or F^O mole fraction has been 
introduced, along with x\ - 1 - h. 

Simple though it appears, Eq. (11) has sev- 
eral interesting and critical features that are 
worth highlighting. In the first place, only V-T 
self-deexcitation of N2 and the V-VR exchange 
reaction of Eq. (3c) survive the cancellation of 
terms caused by application of Eq. (10). 
Secondly, since k 3 is about six orders of 
magnitude larger than k\ (at least around room 
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temperature according to acoustical absorption 
data), very little water vapor is required to mask 
the k\ contribution. Accordingly, it is not 
difficult to imagine the relaxation in most 
laboratory systems being controlled entirely by 
V-VR collisions with H 2 O molecules. A third 

feature is the dependence of pi \ on as well as 
T, which is never possible for single-species 
LHO’s and does not occur here unless 
temperatures are sufficiently high that second- 
order population terms are important (since the 

e* 2 multiplies the energy difference in Eq. (8)). 

A more convenient representation is 
obtained by assuming the relatively 

instantaneous establishment of QSS Boltzmann 
distributions for the populations of vibrational 
levels of each species. This is achieved by fast, 
near-resonant V-V exchange reactions within 
species and corresponds to relaxation 

proceeding through a sequence of such 
distributions. These QSS's are also the basis for 
defining species vibrational temperatures 
according to 



with the translational temperature T replacing 
T v j for the equilibrium energies e- . 

Unless the moist N 2 gas is extremely dilute, 
the existence of a vibrational temperature for N 2 
seems practically assured, as does the existence 
of a vibrational temperature for H 2 O because of 
the very efficient V-V transfer collisions between 
strong dipoles effectively compensating for the 
relative rarity of H 2 O-H 2 O collision events. On 
the other hand, the possibility that the two vibra- 
tional temperatures are equal to each other is 
somewhat remote because the interspecies V-VR 
collisions in reaction (3c) are not nearly as 
probable as the intraspecies V-V collisions 
involved here. The result is that current 2T CFD 
codes (7) T VI = T v j for all i and j) may be less 
adequate than multiple-T codes in which each 
vibrational temperature is allowed to be different 
or perhaps not to exist at all. 

In any event, for the very special gas treated 
in this paper and because of the aforementioned 

QSS for e* 2 culminating in Eq. (7), we can use 


Eq. (7) again to express e* 2 in terms of parame- 
ters mostly pertaining to N 2 . This procedure 
reintroduces some of the V-T reactions that can- 
celed when deriving Eq. (11). The final results 
for small tbO mole fraction li can be written as 
follows, with T v = T V/i = T v ( N 2 ), t>i = A(N 2 ), 
and f >2 = (HtbO): 



where the important correction factor for 
moisture content satisfies 



Expanding versus Contracting Flows 

Clearly the dominant feature of the moisture 
correction factor in Eq. (16) is the factor hk^lk j 
Although the ratio k^/k] may diminish somewhat 

at high temperatures from its value of about 10 6 
at room temperature, as determined from 
acoustical absorption data, 7 the message seems 
indisputable: traces of water vapor have the 
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potential of dramatically accelerating the 
vibrational relaxation of otherwise pure 
nonequilibrium diatomic gases. Moreover, we 
easily understand this behavior in terms of the 
extremely efficient, near-resonant V-VR energy 
transfer collisions between N2 molecules and 
tbO molecules, as opposed to the very 
infrequent V-T collisions of the diatomic 
molecules amongst themselves. We also know 
that the relatively large £3 results from the 
unique geometry of tbO molecules that 
produces a near-continuum of rotational levels 
superimposed on each vibrational level to enable 
resonances with almost any excited diatom. 

The remaining difficulty, of course, is that 
the potentially very large correction factor 
applies to contracting as well as to expanding 
flows. Why is it, then, that recent literature 
dwells almost entirely on anomalously fast 
relaxation in the latter? A possible answer is that 
this has not always been so. Even though the 
physics underlying the effects of water vapor 
was not completely understood, the importance 
of this impurity was recognized empirically and 
early-on by researchers attempting to extract 
relaxation times from shock-tube data. In fact, 
the high-temperature relaxation data base that 
we use today is almost exclusively that derived in 
the 1960's by researchers who had learned to 
take great care in purifying their gases. It is 
obvious, therefore, that huge discrepancies can 
result if these same relaxation times are used in 
comparing theoretical predictions for pure gases 
with data from expanding flows in reflected 
shock tunnels in which similar care is not taken 
in removing impurities. Such theory will always 
underestimate the relaxation rates observed in 
expanding flows. 


Concluding Remarks 

We have developed the rationale and derived 
simple expressions to show the remarkable 
effects of trace amounts of water vapor in 
increasing vibrational relaxation rates in 
nonequilibrium diatomic gases. Since such 
increases arc predicted for both expanding and 
contracting flows, we have further proposed that 
the anomalous relaxation observed in the former 
may be caused by greater water content in 
reflected shock tunnels than in the shock tubes 
from which most relaxation times were deduced. 
This explanation is in marked contrast to the 
more conventional one based on the effects of 
anharmonicity of real molecules. 
Anharmonicity effects do distinguish between 


expansion and contraction, but fall far short of 
accounting quantitatively for the differences 
observed. 

In any event, a case is made for the impor- 
tance of including impurities like tbO in 
modeling flow fields in laboratory test facilities 
in which nonequilibrium properties and the 
internal states of constituent molecules play 
essential roles. Examples range from possible 
effects of molecular excitation on combustion 
chemistry in scramjet engines to the availability 
of excitation energy for conversion into flow 
energy in supersonic and hypersonic nozzles. 
Applications of the theoretical framework 
provided depend, of course, on the availability 
of an adequate set of rate coefficients as 
functions of temperature. These are in short 
supply, especially for V-V and V-VR collisions. 
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